Ground-state properties of even-even N = Z nuclei within the Hartree— Fock— BCS and 

Higher Tamm— Dancoff approaches 
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We calculate the ground-state properties of well deformed, even-even N — Z nuclei in the region 
between 56 Ni and 100 Sn within two different approaches, focusing on the binding energy and defor- 
mation and pairing properties. First, we employ the Hartree-Fock-BCS (HFBCS) approximation 
with the Skyrme effective nucleon-nucleon interaction and discuss how the results depend on the 
parameterization of the interaction and on the pairing force parameters adjusted in various schemes 
to reproduce the experimental odd-even mass differences. Then, within the Higher Tamm-Dancoff 
Approximation (HTDA), which explicitly conserves the particle number, we calculate the same 
properties starting from the HFBCS solutions. The HTDA treatment of the ground-state correla- 
tions is converged within a n-particle-n-hole expansion using up to n=4 particle-hole excitations of 
the pair type (in the sense of Cooper pairs). We compare the ground-state properties calculated 
in these two descriptions of pairing correlations and deduce the importance of the particle-number 
conservation in weak pairing regimes. Finally, we extend the HTDA calculations so as to include 
the proton-neutron residual interaction and investigate the role of proton-neutron pairing on the 
above ground-state properties. 

PACS numbers: 21.60. Jz,21.60.Cs 



I. INTRODUCTION 

Nuclei with an equal number of neutrons and protons 
are of a special interest in many respects. The similarity 
of the neutron and proton single-particle (sp) spin-space 
wave functions, in the vicinity of the chemical potentials, 
allows for rather interesting physical studies associated 
with isospin mixing, fi decays and proton-neutron cor- 
relations. Such peculiar correlations are encountered in 
pairing properties either in the T — or T — 1 chan- 
nels. Another intriguing property is the so-called Wigner 
term in a liquid-drop approaches, corresponding to a very 
sharp pattern of extra stability close to N = Z. These 
nuclei are also of importance for the rapid-proton astro- 
physical process. Up to A ~ 60 for even values of N = Z, 
these nuclei are very close indeed to the proton drip line. 
However, contrary to the situation near the neutron drip 
line and due to the Coulomb barrier, their description in 
a mean-field-plus-correlations approach is not marred by 
the necessity of dealing accurately with the treatment of 
the continuum. If one is interested in specific properties 
of such nuclides and their neighbors (such as the Wigner 
term) , it is desirable to use a description of pairing corre- 
lations which is not blurred by fluctuations of the neutron 
and proton numbers as is the case when using Bogoliubov 
quasiparticle vacua as ansatz for the ground-state (GS) 
wave functions. 

In this paper we are concerned with the mean-field and 
beyond-mean-field descriptions of well-deformed even- 
even N = Z nuclei lying between the doubly magic nuclei 
56 Ni and 100 Sn. To do so we use various Skyrme effec- 



tive interactions (SIII [l[ and SLy4 [2[) supplemented by 
two different treatments of the pairing correlations: a 
la Bogoliubov- Valatin (BCS wave function) and within 
the particle-number conserving approach dubbed as the 
Higher Tamm-Dancoff Approximation (HTDA) Q, in 
both cases for like-particle pairing correlations. Then we 
add in the latter treatment proton-neutron pairing corre- 
lations in T = or T = 1 isospin channels. 

In view of its widespread use for several decades, the 
BCS approach for like-particle pairing properties is not 
reviewed here. Only some practical details (among which 
the choice of the relevant average matrix elements is of 
paramount importance) are discussed. In contrast the 
HTDA approach is less widely known and is briefly de- 
scribed here. 

The HTDA approach may be presented as a treatment 
of correlations in a highly truncated shell model whose 
practicability and efficiency rely on the fast convergence 
of the particle-hole expansion. This has been shown to 
be realized upon choosing for the particle-hole quasipar- 
ticle vacuum a relevant Hartree-Fock (HF) solution as- 
sociated self-consistently with the one-body reduced den- 
sity matrix of the correlated wave function. The HTDA 
approach was applied for the first time to describe the 
ground and isomeric states in 178 Hf Q and then odd nu- 
clei and most general isomeric states Q. A Routhian- 
HTDA scheme was then proposed to describe the su- 
perdeformed yrast bands in the A ~ 190 region A 
preliminary HTDA study of the GS pairing correlations 
in 64 Ge, including isovector and isoscalar residual interac- 
tions, was presented in Ref. @. Here the HTDA approach 
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is applied for the first time in systematic calculations of 
the properties of medium-mass, proton-rich nuclei. 

The outline of the present work is the following. In 
Sect. [XT] we describe the theoretical background (mostly 
the HTDA formalism) and explain at length in Sec. IIIII 
how the calculations are carried out in practice, espe- 
cially regarding the pairing strength fitting strategy in 
both BCS and HTDA approaches and the optimization 
of the harmonic-oscillator basis parameters in the HFBCS 
calculations. The presentation and discussion of the ob- 
tained results is organized in three steps. First, in Sec. II VI 
we focus on the GS properties obtained within the HF- 
BCS approach, and in Sect. [V] on the results obtained 
within the HTDA formalism without proton-neutron cou- 
pling in the particle-particle channel. Then we carry out 
in Sect. I VII a comparative study of pairing properties de- 
pending on the pairing treatment and the fitting scheme. 
Finally we extend in Sec. I VIII the HTDA calculations so 
as to include the full isovector and isoscalar pairing in- 
teraction. The main conclusions of our study are drawn 
in Sec. IVTHl 



II. THEORETICAL BACKGROUND 

In the HFBCS approach we use the Skyrme effective 
interaction in the particle-hole channel and the senior- 
ity force in the particle-particle channel, and expand the 
single-particle wave functions in the cylindrical harmonic- 
oscillator basis, as detailed in Ref. [3]. 

We focus in the remainder of this section on the HTDA 
approach. Its purpose is to describe various nucleon cor- 
relations (such as pairing and RPA) on the same foot- 
ing. Among the many particle-hole excitations on a 
Slater determinant vacuum (noted here l^o) 1 ); pair ex- 
citations around the Fermi surface play an essential role. 
By construction the HTDA approach is an extension 
of the Tamm-Dancoff approximation to higher order of 
particle- hole excitations, so it may be regarded as a trun- 
cated shell model. The rapidity of the convergence of 
the particle-hole expansion, and thus the tractability of 
the approach, depends on the realistic character of 
A fast convergence is expected to be reached when the 
quasi-vacuum is defined self-consistently in such a way 
that some many-body effects of the correlations are al- 
ready taken into account at the mean-field level. 

Let us consider the effective Hamiltonian 

H = K + V, (1) 

where K denotes the kinetic energy operator and V an ef- 
fective interaction. For the wave function |$o) we choose 



the HF solution, i.e., the eigenstate of the Hamiltonian 
Hup defined by 

H HF = K + V>hf (2) 

where the potential Vjjf denotes the one-body reduction 
of V for |<E>o) and is self-consistently obtained from the 
many-body reduced density matrix p of the correlated 
solution Y$>) for the desired number of particles. We have 
thus 

#hf|*o) = #o|*o), ( 3 ) 

where Eq is the associated eigenenergy and 

Hn F = K+ Vhf • (4) 

This approach also allows to include various one-body 
constraints (on the nuclear deformation for example) in a 
simple way since the constraint operator can be absorbed 
in the definition of -Hhf- I n addition we assume here 
that the GS solutions possess the time-reversal, axial and 
parity symmetries. 

The quasi-vacuum |<5> ) may now serve to construct 
an orthonormal many-body basis in which we diagonal- 
ize the Hamiltonian H. In principle to build this ba- 
sis we should include, besides |$o) = \&o) ® 1*^0 )> tnc 
particle-hole excitations of all orders (from 1 to the par- 
ticle number) created on |$o), noted generically |$ n ) for 
n-particle-n-hole (np-nh) excitations. The total GS wave 
function \^f) = \^ T ) <8> \ ) can therefore be decomposed 
in the following way 

|*)=Xoo|*o)®l*o')+ 2 Xio|*D®|*o') 
(i P -ih) T 

(lp-lh) T , (lp-lh) T 

(lp-lh) T , 

(2p-2h) T (2p-2h) T , 

H 

(5) 

where r and t' denote two different charge states. How- 
ever practical calculations require to truncate this expan- 
sion. Based on the former studies in the HTDA frame- 
work [3, |4j, |6( , we may assume that the components of the 
pair-excitation type dominate in the GS solution. The 
set of products of Slater determinants |<E>[) <g) |$J ) is an 
orthonormal basis of the physical space accessible to a 
nucleus having N neutrons and Z protons. Assuming 
time-reversal symmetry, the coefficients Xi m Eq. (0) are 
real and, if we take |^) normalized to unity, obey the 
relation 



1 For the sake of clarity in the notation, we reserve the letter <J> for 
a Slater determinant and the letter for a correlated state. 
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(6) 
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It can be easily shown that the expression {5} of en- 
sures that |\&) is an eigenstate of the particle- number op- 
erator: N\ty) = with A = N+Z. Finally, to obtain 
the correlated ground state we diagonalize the Hamilto- 
nian defined in Eq. ([1]) in the retained many-body basis. 
Now, let us rewrite the Hamiltonian H as 

H = ($ \H\$ )+H lQF + V Ies , (7) 

where the independent quasiparticle Hamiltonian -f/iQp 
and the residual interaction V^ cs are defined by 

i?IQP=-ff H F-($o|^HF|$o), (8) 

V Ies = V-VHF + {^o\V\^o) ■ (9) 

It is interesting to note that these definitions give a van- 
ishing expectation value of Hiqp and V les for the HF so- 
lution |$o). The independent quasiparticle Hamiltonian 
can also be expressed as 

Hiqp = ^2^kVkVk (10) 

k 

where is equal to the energy e^. of the single-particle 
state |fe) if \k) is a particle state with respect to |<&q) 
or equal to — ej, if \k) is a hole state. In Eq. (fTQ|) . ryl is 
the creation operator at associated with \k) if \k) is a 
particle state or the annihilation operator at if \k) is a 
hole state. The matrix element of H between two Slater 
determinants |$,) and of the multi-particle multi- 
hole basis therefore takes the form 

(QilH^j) = %«$o|^HF|$o)+^4i T) ) + < $ il^es|^), 

T 

(11) 

where is the particle-hole excitation energy associ- 

ated with \&J) and calculated with respect to the vacuum 
1*5) as 

e% ] = E ■ ( 12 ) 

where the sum runs over all the single-particle states \k) 
contained in \<f>J). Since the residual interaction con- 
tains only one- and two-body operators, the matrix ele- 
ment ($j|V^es|^j) vanishes when and |$A differ by 
particle-hole excitations of order three or higher. It can 
be calculated in terms of two-body matrix elements of V 
by application of the Wick theorem generalized to quasi- 
vacua as shown in Ref. [ij. 

III. CALCULATION PROCEDURE 

Since we are in a region of nuclear deformation instabil- 
ity (shape coexistence) our final results can be sensitive 



to the choice of the effective interaction and to the pair- 
ing treatment. This is why two parameterizations of the 
phenomenological Skyrme interaction, namely the SIII 
and SLy4 ones are used. In the particle-particle channel 
a simple seniority ansatz is adopted in the BCS calcula- 
tions. It is specified by the value G of the constant pairing 
matrix elements between any single-particle states in the 
canonical basis retained in this part of the calculations. 
The value of G is adjusted to reproduce the so-called em- 
pirical pairing gaps evaluated in the 3-point or a 5-point 
formula Q. The differences between these two experi- 
mentally deduced quantities are large in this region and 
so are the experimental errors for nuclear masses. In ad- 
dition the meaning of the pairing gap derived from finite 
mass differences is no longer clear in the N = Z case. 
Therefore, we find it necessary to make a comparison 
of the results obtained in the two cases which represent 
stronger (5-point fit) and somewhat weaker (3-point fit) 
pairing regimes. We describe in the next subsection the 
fitting procedure for G as well as the harmonic-oscillator 
(HO) basis parameters optimization carried out in the 
HFBCS framework. The obtained basis parameters are 
then used in the HTDA calculations assuming that they 
do not differ significantly from the values that would re- 
sult from an optimization in the HTDA approach. 

Once the solutions corresponding to the equilibrium 
deformations determined in the HFBCS calculations are 
found, we perform perturbative HTDA calculations (one 
diagonalization of the HTDA matrix is performed on top 
of the HFBCS calculations) without any constraints on 
the deformation. Since we are interested in the GS corre- 
lations of even-even nuclei, the many-body basis includes 
here only pair excitations, i.e., excitations where nucle- 
ons occupying twofold Kramer-degenerate hole levels are 
scattered into Kramcr-dcgcncratc particle levels. As for 
the residual interaction in the HTDA approach, we choose 
a S force whose strength is adjusted as explained in sub- 
section Bl In this way our formalism not only retains 
the realistic character of self-consistent mean-field calcu- 
lations, but also ensures the particle number conserva- 
tion. 



A. Pairing strength adjustment in the HFBCS 
approach 

We shall follow here the steps of Bonche and collab- 
orators [9j and apply the BCS approximation with the 
seniority pairing interaction. For the sake of defining our 
notation and for completeness, we recall the relevant ex- 
pressions involved, omitting the isospin r index. First, 
the seniority antisymmetrized matrix element is given by 

Vkik'l' = ~G fkfk'S a 6 llW (13) 
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with the smooth cut-off factor 



fk — 



1 + cxp[(e fc - e F - Ae)//j] 



(14) 



Here as well as in the HTDA case, the Fermi level ep is 
defined by 



1 , . 



(15) 



where e n and e„+i denote the energies of the last occupied 
and the first empty single-particle state, respectively (in a 
pure HF picture), and Ae and \x denote the cut-off energy 
and the diffuseness parameters. The pairing gap can be 
expressed as 

A fc = f k A (16) 
where the state independent gap A is given by 



A= G y /|A 



k>0 Eqp 



with the quasiparticle energy E^) defined by 

4$ = >J% + fZ* 2 - (is) 

In Eq. p7p, the sum runs over all the pairs of Kramer 
degenerate single-particle states, but in fact the cut-off 
factor f}. suppresses the contributions of single-particle 
states lying at least about Ae + \x above the Fermi level. 
From the average particle number conservation 



(17) 



fc>0 



(19) 



we can deduce the expression of the chemical potential A 



N ■ 



k>0 \ qp 



(20) 



Finally the pairing energy takes here the simple form 



F - — 

-C'pair — q 



(21) 



We now discuss the determination of the pairing 
strengths G^ related to the actual matrix elements G^ 
through the following prescription 



G (t 



G, 



11 + A^r 



(22) 



First we approximately take into account the Coulomb 
reduction of proton pairing by assuming that, as the 



Hartree-Fock-Bogoliubov (HFB) calculations with the 
Gogny D1S interaction tend to indicate [To( | : 



G 



(p) 



0.9 G 



(n) 



(23) 



Then, for a given parameterization of the Skyrme force 
and a given set of basis parameters {No , b and q in the no- 
tation of Ref. p}), we determine the GS deformation fa of 
each nucleus using a reasonable initial G- value, assumed 
to be the same for all the nuclei under study. Then we de- 

(n) 

duce G — Gq from a least-square fit to the experimental 
minimal quasiparticle energies through the 3-point and 5- 
point formulae [j| [ll[ and using the same single-particle 
spectrum (the one for the charge state r corresponding to 
the converged solution at fa)- We thus have to minimize 
the following function 



X 2 (G) 



1 



2iVn 



JVnucl 

E E 

i—l r—n,p 



[£r(G)].-[^l 



where £ T denotes the lowest quasiparticle energy of the 
nucleons of type r and i is an index running over the A^ nuc i 
nuclei included in the fit. With the obtained G- value, we 
determine the new GS deformation of each nucleus and 
minimize again ^{G) to find a new G- value. This pro- 
cedure is repeated until the simultaneous convergence of 
G and the GS deformations is reached. In practice, 
it is necessary to scan a wide range of deformations to 
find the lowest local minimum of the deformation energy 
curve. The latter is determined using the basis parame- 
ters b = \J ttiluo / h, with luq = (cj^Cl^) 1 / 3 , and q = lq±/uj z 
deduced from an approximate expression for a HO po- 
tential (see Ref. Q). This approximate optimization re- 
quires only the knowledge of bo, the optimized value of 
b{q) for a spherical solution. We carry out this study with 
b a = 0.505, which is approximately the actual optimal 
value for the considered nuclei. By varying the basis pa- 
rameters in the calculated ground states of all nuclei, we 
have checked that the optimal G-value does not change 
significantly. The obtained values of G are reported in 
Table U 



B. Pairing strength adjustment in the HTDA 
approach 

As mentioned earlier, in the HTDA calculations we use 
a S force in the particle-particle channel. Both the cou- 
pling constant and the cut-off parameter are necessary to 
define fully the interaction. In our approach it is then 
necessary to fix the strength of this interaction and this 
is done by adjusting it so as to reproduce physical quanti- 
ties for the considered nuclei, e.g., the phenomenological 
gaps. For that purpose, assuming that the appearance of 
the pairing gap is related to the breaking of the Cooper 
pair of lowest energy, we block in the HTDA calculations 
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TABLE I: Empirical pairing gaps deduced from the odd- 
even mass differences (using the Atomic Mass Evaluation 
AME2003 H2|) in 5-point and 3-point formulae [EEH, except 
for A<, 5) ( 80 Zr) for which we use the nuclear mass calculated 
by Moller and Nix [l3[|. The standard estimate of the pairing 
gap is given for comparison. In the last column the 

experimental binding energy per nucleon is also indicated. All 
values are in MeV. 



Nucleus Ai 5) A£ 5) Ai 3) A ( p 3) 12/ y/A E/A 



64 Ge 


2.07 


1.85 


1.48 


1.14 


1.50 


8.5294 


68 Se 


2.24 


2.01 


1.66 


1.37 


1.45 


8.4773 


72 Kr 


1.72 


1.73 


1.24 


1.10 


1.41 


8.4293 


76 Sr 


1.46 


1.635 


0.88 


1.07 


1.37 


8.3938 


80 Zr 


1.94 


1.57 


1.31 


0.99 


1.34 


8.3741 



the single-particle level (neutron or proton) closest to the 
Fermi energy. Then, we adopt the difference between the 
expectation value of V^ es in a normal (n) and blocked (b) 
calculations as a proper measure of the pairing correla- 
tions that can be compared to the experimental odd-even 
mass differences. Namely, we define 

A = - fg.) - (E^ Eg},) , (25) 
where E = and £t QP = (*|.Hiqp|*). 

IV. RESULTS OF THE HARTREE— FOCK— BCS 
CALCULATIONS 

The five selected well-deformed N = Z nuclei are 64 Ge, 
68 Se, 72 Kr, 76 Sr and 80 Zr. For one of the two consid- 
ered effective interactions (SIII) we use several pairing 
windows with different values of Ae (6, 8 and 10 MeV) 
and n (0.2 and 0.5 MeV). In the fitting process, whereas 
the choice of Ae is rather unimportant in the considered 
range, we find that it is not the case for fi. For instance, 
with ll = 0.5 MeV and Ae = 6 MeV, the iterative proce- 
dure to adjust of G undergoes oscillations preventing to 
reach convergence, contrary to all the other pairing win- 
dows considered here. We thus retain for further calcula- 
tions the values Ae = 6MeV and /i = 0.2MeV stemming 
as the best choice from the fit to the 5-point experimental 
gaps. The optimal pairing strengths G opt for the different 
Skyrme interactions and fitting schemes are displayed in 
Table |TT] together with the root-mean-square error on the 
experimental quasiparticle energies ex a 

<ta = yVtGopt) ■ (26) 



TABLE II: Optimal G-values and corresponding root-mean- 
square error on the experimental gaps cta (in MeV) obtained 
with the SIII and SLy4 Skyrme interactions using the 3-point 
and the 5-point formulae with Ae = 6 MeV, /i = 0.2 MeV and 
No = 10. 



Force Formula G op t 


CTA 


SIII 


3-point 


17.7 


0.120 


SIII 


5-point 


20.6 


0.216 


SLy4 


3-point 


17.2 


0.194 


SLy4 


5-point 


19.9 


0.260 



These adjustments have been performed with a basis 
size defined by N = 10. The same adjustment proce- 
dure has also been carried out with a much larger HO 
basis (Nq — 16) in the illustrative case of the pairing 
window parameters Ae = 6 MeV and \x = 0.2 MeV. The 
optimal value of G and the associated root-mean-square 
error turns out to be very close to that obtained with 
N Q = 10. This justifies the choice of N = 10. 

It is worth adding that the authors of Ref. Q per- 
formed similar calculations in the same mass region (in- 
cluding for 76 Sr and 80 Zr) with the SIII Skyrme interac- 
tion and the following set of parameters for the seniority 

force: Ae = 5 MeV, fi = 0.5 MeV, G { Q n) = 13.5 MeV and 

(v) 

GJp = 16.5 MeV. The pairing strengths were determined 
from the experimental quasiparticle energies extracted in 
the same way as we did from the experimental binding 
energies. 

With the optimal values of the pairing strength of Ta- 
blc|TT]we calculate several GS properties related to the nu- 
clear size (through the root-mean-square mass radius r m ) 
and deformation (through 02 and the mass quadrupole 
Q20 and hexadecapole Q40 moments), the binding en- 
ergy per nucleon (E/A) as well as pairing quantities (the 
BCS pairing gaps A„ and A p , the chemical potentials X n 
and A p , and the minimal quasiparticle energies £ n and 
£ p ). The definitions of P2, r m , Q20 and Q40 can be found 
in the Appendix. The obtained results are reported in 
Table Mil The most striking difference between the two 
Skyrme interactions is that they yield very different GS 
deformations for 80 Zr (strongly prolate with SIII, spher- 
ical with SLy4). Moreover the GS deformation of 76 Sr 
drastically depends on the pairing strength when using 
the SLy4 interaction. 

Since the GS deformations calculated with the SIII in- 
teraction are in agreement with the experimental data, 
especially for 76 Sr Q and 80 Zr p]|, we perform HTDA 
calculations only with this interaction. 
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TABLE III: Ground-state properties obtained with the SIII and SLy4 Skyrme interactions and the pairing strengths obtained 
through the 3-point and the 5-point formulae adjustment procedure. All the quantities in the columns at the right of Q40 are 
expressed in MeV. 



Force Formula Nucleus (3 2 r m (fm) Q20O) Q4o(b ) E/A A„ A p 



An 



At, 



SIII 
SIII 
SIII 
SIII 
SIII 
SIII 
SIII 
SIII 
SIII 
SIII 



3-point 
3-point 
3-point 
3-point 
3-point 



64 Ge 
68 Se 
72 Kr 
76 Sr 
80 Zr 



0.200 
-0.267 
-0.340 
0.390 
0.401 



3.917 
4.015 
4.117 
4.238 
4.318 



2.651 
-3.324 
-4.550 
7.498 
8.441 



0.0097 
0.0535 
0.1038 
0.2353 
0.1879 



8.4216 
8.3698 
8.3229 
8.3024 
8.2592 



1.375 
1.173 
0.966 
0.425 
0.967 



1.138 
0.891 
0.579 
0.000 
0.623 



12.642 
13.040 
13.057 
13.653 
13.131 



2.713 
-2.760 
-2.479 
-2.584 
-1.586 



1.376 
1.502 
1.299 
1.132 
1.185 



1.140 
1.222 
1.081 
0.971 
0.965 



5-point 
5-point 
5-point 
5-point 
5-point 



64 Ge 
68 Se 
72 Kr 
76 Sr 
80 Zr 



0.191 
-0.244 
-0.259 
0.381 
0.384 



3.918 
4.012 
4.089 
4.234 
4.308 



2.515 
-3.062 
-3.546 
7.301 
8.031 



0.0141 
0.0460 
0.0544 
0.2133 
0.1995 



8.4485 
8.3967 
8.3472 
8.3087 
8.2741 



2.035 
2.074 
2.064 
1.380 
1.766 



1.690 
1.709 
1.710 
0.947 
1.387 



12.775 
13.119 
13.228 
13.537 
13.256 



-2.830 
-2.807 
-2.477 
-2.515 
-1.685 



2.036 
2.238 
2.065 
1.683 
1.828 



1.690 
1.876 
1.712 
1.330 
1.495 



SLy4 

SLy4 

SLy4 

SLy4 

SLy4 

SLy4' 

SLy4 

SLy4 

SLy4 

SLy4 



3-point 
3-point 
3-point 
3-point 
3-point 



64 Ge 
68 Se 
Kr 



72 



71, 



Sr 
Zr 



0.166 
-0.256 
-0.169 
0.392 
0.000 



3.887 
3.996 
4.047 
4.218 
4.158 



2.141 
-3.172 
-2.359 
7.470 
0.000 



0.0010 
0.0344 
-0.0038 
0.2380 
0.0000 



8.4719 
8.4347 
8.3741 
8.3402 
8.3377 



1.162 
0.672 
1.275 
0.000 
0.461 



0.978 
0.133 
1.058 
0.000 
0.000 



12.687 
12.984 
13.237 
13.627 
13.911 



-2.831 
-2.768 
-2.541 
-2.597 
-2.352 



1.162 
1.551 
1.295 
1.226 
1.458 



0.979 
1.322 
1.075 
1.149 
1.238 



5-point 
5-point 
5-point 
5-point 
5-point 



64 Ge 
68 Se 
72 Kr 
76 Sr 
80 Zr 



0.150 
-0.228 
-0.164 
0.000 
0.000 



3.887 
3.992 
4.048 
4.097 
4.158 



1.922 
-2.844 
-2.297 
0.001 
0.000 



0.0035 
0.0288 
-0.0010 
0.0000 
0.000 



8.4908 
8.4489 
8.3958 
8.3614 
8.3450 



1.687 
1.736 
1.839 
1.805 
1.365 



1.426 
1.376 
1.533 
1.530 
1.065 



12.764 
13.064 
13.296 
13.952 
13.931 



2.907 
2.823 
2.608 
2.863 
2.383 



1.693 
2.153 
1.852 
1.867 
1.971 



1.430 
1.837 
1.543 
1.572 
1.655 



V. RESULTS OF HTDA CALCULATIONS 

A. Convergence of the HTDA solutions 

In the HTDA calculations two types of truncation 
schemes need to be defined. The first one is concerned 
with the maximal order in the many-particle many-hole 
basis, the second one with the single-particle states from 
which this basis is built. In the latter case, we are fac- 
ing thus a situation met in customary BCS calculations. 
Typically, we limit our single-particle subspace to the 
configuration-space window defined as 

|e fc - ei r|<Ae, (27) 

where ep is the Fermi energy and Ae is a cut-off energy. 
The actual value of Ae is chosen such that single-particle 
states left out would not contribute significantly except 
by a renormalization of quantities measuring the inten- 
sity of pairing correlations (such as the correlation energy 
defined below). As already used in HTDA calculations 
for the 64 Ge nucleus Q, we retain here Ae = 12 MeV 
for both charge states. As in the BCS treatment, the 
two-body matrix elements arc multiplied by the smooth 
cut-off factor of Eq. (J33) with Ae = 12 MeV and fj, = 0.2. 



A detailed study of the particle-hole expansion of the 
HTDA ground state was performed in Ref. [l6| in the 
picket-fence model with 8 or 16 levels filled with 8 or 16 
particles, respectively. The authors obtained the conver- 
gence of GS solutions toward the exact solution of the 
Richardson model for a 6p-6h space in the standard nu- 
clear pairing regime. Nevertheless, it may be expected 
that the realistic character of the HF vacuum calculated 
with a realistic effective interaction ensures a faster con- 
vergence of the particle-hole expansion. This is what the 
studies of Refs. using a 4p-4h space tend to indicate. 

Within the above truncation scheme for the 
configuration-space window obtained in the HFBCS 
calculations with the SIII force and G = 20.6 MeV, 
we study the convergence of the HTDA solutions as a 
function of the many-body basis size. The three quan- 
tities under consideration are i) the correlation energy 
defined as the difference between the expectation values 
of the Hamilton operator evaluated in the correlated and 
uncorrelated states 

E COII = {V\H\V) - (4 |^|*o) , (28) 

ii) the occupation probabilities vf which are defined in 
the HTDA approach as the diagonal matrix element of 
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FIG. 1: Convergence of the total correlation energy for Kr 
as a function of the pairing strength Vo. The percentage dif- 
ference between correlation energies obtained in many-body 
spaces that differ by 2p-2h are shown. 

the one-body density p in the single-particle basis 

v* = p ii = ma\a i \*), (29) 

where a\ and are the creation and annihilation op- 
erators associated with the single-particle state \i), re- 
spectively, and iii) the mass quadrupole moment (see the 
Appendix for definitions). Since the one-body density 
is not diagonal in the HF basis a transformation to the 
canonical basis is done to obtain the vf values. 

The percentage difference of E con . between solutions 
obtained in 2p-2h, 4p-4h and 6p-6h spaces of pair excita- 
tions are shown for 72 Kr as a function of the strength Vo 
of the residual S interaction in Fig. [TJ The difference be- 
tween 2p-2h and 4p-4h solutions is found to be very large 
and increases nearly linearly with the pairing strength 
Vq. The discrepancy between 4p-4h and 6p-6h solutions 
reaches 5% in the strong pairing region. A similar behav- 
ior is found for the other nuclei. 

The occupation probabilities of single-particle levels 
obtained in the different spaces and with different Vq- 
values are shown in Fig. [2] for 76 Sr. There are no con- 
spicuous differences between the various solutions except 
in the stronger pairing case, where the difference between 
the vf values in 2p-2h and larger spaces becomes substan- 
tial. Therefore, the calculated quantities like quadrupole 
moments and radii converge quickly with the maximal 
order of particle- hole excitations. 

In Tables IIVI and [V] the values of correlation energies 
and quadrupole moments obtained in 2p-2h, 4p-4h and 
6p-6h spaces are given for all nuclei and for two real- 
istic values of the strength of the residual interaction, 
namely V = 400 MeV.fm 3 and V = 320 MeV.fm 3 (see 
the next subsection for details). The largest discrep- 
ancy between 4p-4h and 6p-6h results occurs for 72 Kr 
with Vq = 400 MeV.fm 3 where the difference for the cor- 
relation energy reaches 3.6%. With the pairing strength 



TABLE IV: Correlation energy and mass quadrupole moments 
obtained in the calculations using 2p-2h, 4p-4h and 6p-6h 
spaces with the 5-force strength Vo = 400 MeV.fm 3 . 



Nucleus 


£ corr (MeV) 
2p-2h 4p-4h 6p-6h 


2p-2h 


Q 2 o(b) 
4p-4h 


6p-6h 


64 Ge 


-4.763 


-6.383 -6.556 


2.670 


2.645 


2.641 


68 Se 


-4.127 


-5.060 -5.205 


-3.301 


-3.306 


-3.306 


72 Kr 


-4.306 


-5.348 -5.594 


-4.575 


-4.528 


-4.519 


76 Sr 


-4.010 


-4.996 -5.042 


7.493 


7.482 


7.480 


80 Zr 


-4.077 


-5.046 -5.228 


8.453 


8.449 


8.448 



TABLE V: Same as in Table [TV] with the strength V = 
320 MeV.fm 3 of the 5-force. 



Nucleus 


Scorr(MeV) 

2p-2h 4p-4h 6p-6h 


2p-2h 


Q 20 (b) 
4p-4h 


6p-6h 


64 Ge 


-2.992 


-3.796 -3.870 


2.676 


2.656 


2.655 


68 Se 


-2.404 


-2.798 -2.842 


-3.301 


-3.304 


-3.304 


72 Kr 


-2.558 


-2.936 -3.026 


-4.595 


-4.574 


-4.571 


76 Sr 


-2.493 


-2.696 -2.758 


7.498 


7.493 


7.493 


80 Zr 


-2.421 


-2.768 -2.825 


8.456 


8.456 


8.456 



Vo = 320 MeV.fm 3 the discrepancies for _E C orr amount 
on the average to 1%. For all five nuclei the differences 
between the values of the quadrupole moments obtained 
in 4p-4h and 6p-6h calculations are negligible. 

These results suggest that it is reasonable to truncate 
the particle-hole expansion at order 4. As a matter of 
fact, this allows for a satisfactory accuracy of the calcu- 
lations at a rather low cost in terms of computation time. 
A typical number of many-body configurations to be han- 
dled for each charge state is about 5000 while it may rise 
typically up to 50000 or higher when including 6p-6h exci- 
tations. It is worth adding here that in these calculations 
no additional cut-off on the maximal particle-hole energy 
of the scattered pairs is applied. The addition of such a 
cut-off for the particle-hole excitation energy would yield 
an even faster convergence without losing any substantial 
accuracy (l6j and may be one way to handle calculations 
of higher particle-hole excitation order. 



B. Adjustment of the strength of the residual 
5-interaction 

The fit of the strength of the S interaction in the HTDA 
approach is performed in the space of up to two-pair ex- 
citations. Two adjustment schemes are considered: one 
where the same strength is retained for neutrons (Vo n ) 
and protons (Vo p ), the other one where Vo p is reduced 
by 10% with respect to Vo n to account for the Coulomb 
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FIG. 2: Convergence of the occupation probability vj with enlarging of the particle-hole excitations space for neutrons and 
protons in 76 Sr. The cases of three different intensities Vo of pairing interaction are shown. 



anti-pairing effect as in the HFBCS calculations. 

In Table [VTl we present the optimal Vb p -values, noted 
VJ^ pt , and the associated gap root-mean-square deviations 
<ta in each adjustment scheme and using the 3-point and 
5-point formulae. Since the reduction of V^ vt slightly 
improves the quality of the fit, we choose this adjustment 
scheme in further HTDA calculations. 



TABLE VI: Optimal values V opt (in MeV.fm 3 ) and corre- 
sponding root-mean-square errors <ta (in MeV) on pairing 
gaps obtained with the SIII Skyrme interaction using the 3- 
point and the 5-point formulae. 



Formula 


Vb p = 

T/<>Pt 


= Vo„ 


Vbp = 

Tyopt 

v o 


0.9V O n 
CTA 


3-point 


320 


0.278 


340 


0.264 


5-point 


400 


0.285 


420 


0.249 



GS properties are presented in Tables IVTTl and [Villi when 
using the 3-point and 5-point formulae, respectively. 

TABLE VII: Ground-state properties within the HTDA 
framework with the SIII force and the pairing strength ob- 
tained through the 3-point formula adjustment procedure. 
The mass root-mean-square radius r m , quadrupole moment 
and hexadecapole moment are given in fm, barns (b) and b 2 , 
respectively, whereas the binding energy per nucleon E/A and 
the pairing gaps are expressed in MeV. 



Nucleus 


02 


I'm 


Q20 


Qio 


E/A 


A„ 


A p 


64 Ge 


0.201 


3.917 


2.661 


0.0076 


8.4340 


1.801 


1.698 


68 Se 


-0.267 4.014 


-3.315 


0.0530 


8.4004 


1.338 


1.324 


72 Kr 


-0.341 


4.120 


-4.575 


0.1059 


8.3577 


1.164 


1.152 


76 Sr 


0.389 


4.240 


7.494 


0.2344 


8.3344 


2.003 


1.524 


80 Zr 


0.400 


4.320 


8.444 


0.1830 


8.2863 


1.216 


1.011 



C. Ground-state properties obtained within the 
HTDA approach 

As mentioned in Sect. IIV1 the HTDA calculations use 
as a starting point the HFBCS solutions obtained with 
the SIII effective interaction and two sets of values of the 
pairing-strength adjusted to experimental data through 
the 3-point and 5-point formulae. The results for the 



Since we are considering here only the equilibrium de- 
formations as determined by the HFBCS calculations, 
HTDA calculations do not affect much the bulk observ- 
ables, with the exception of the binding energy and re- 
lated quantities such as the pairing gaps. The former 
is indeed much lower (by 2 to 4 MeV) for the HTDA 
solutions than for the HFBCS ones. Regardless of the 
adjustment scheme for the pairing strength, the binding 
energy calculated in the HTDA approach agrees better 
with the experimental values (see Table [XJ) than in the 
HFBCS approach. 
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TABLE VIII: Same as Table IVIIl using the 5-point formula 
adjustment procedure. 



Nucleus 


02 




Q20 


Q40 


E/A 


A„ 


A p 


64 Ge 


0.200 


3.919 


2.651 


0.0089 


8.4723 


2.113 


2.011 


68 Se 


-0.266 


4.017 


-3.318 


0.0536 


8.4326 


2.023 


2.006 


72 Kr 


-0.337 4.121 


-4.532 


0.1027 


8.3898 


1.718 


1.714 


76 Sr 


0.382 


4.242 


7.483 


0.2330 


8.3620 


2.003 


1.524 


80 Zr 


0.400 


4.321 


8.436 


0.1879 


8.3138 


1.513 


1.323 



VI. PAIRING PROPERTIES IN BCS AND 
HTDA APPROACHES 



To evaluate the amount of pairing correlations one 
might think of considering the correlation energy, defined 
as the difference between the expectation values of the 
Hamiltonian in the correlated and uncorrelated (HF vac- 
uum) solutions. However, such a correlation energy has 
no realistic character in practice, in both the HFBCS 
and HTDA approaches. Indeed, there is no consistency 
between the interaction used to generate the mean-field 
and the one building up pairing correlations. Instead, a 
quantity only related to the residual interaction would 
be more significant in this respect. This is the case of 
the so-called condensation energy £" C ond- In the HFBCS 
approach, -E C ond is proportional to the trace of the prod- 
uct of the abnormal density and the pairing field. In the 
HTDA approach, we may define it by 



E, 



cond 



= E, 



- V y 2 E 1 

corr / J Xi ^p- 



h - 



(30) 



where the xf factors are the probability of the configura- 
tion i whose unperturbed particle- hole energy is E ph . 

Another variable that may shed light on the amount 
of pair correlations is the trace of the positive-definite 
operator /3 5 (1 — p)i. Indeed, this quantity expresses the 
non-idempotent character of the density operator p and 
is thus related to correlations. As well known in the BCS 
case, it is related to the abnormal density. Using the 
occupation factor Vi defined in Eq. (|29p . we can write 
the trace of the operator /5 5 (1 — p)* simply as the sum 
J2i UiVi, with Ui = y/l- vf. 

These two measures of pairing correlations calculated 
in the HFBCS and HTDA approaches are compared sep- 
arately for neutrons and protons in Tables IIXI and [Xj 
respectively. The most striking feature is the rather tiny 
variations in the HTDA case from one nucleus to another. 
This is due to the resilience of the HTDA solutions to re- 
act on variations of the level density at the Fermi surface. 
In contrast, it is overemphasized in the HFBCS calcula- 
tions. One example of this is to be found for the protons 
distribution in 76 Sr, where in the HFBCS calculations 



TABLE IX: Neutron condensation energies and diffuseness of 
the neutron Fermi surface obtained in the HFBCS and HTDA 
approaches. The results with SHI force and two fits of the 
strength of the pairing interaction are given in each case. 



Formula Nucleus 



HFBCS 



HTDA 







-E'cond 


Y u i y i 


-E'cond 


Y u i V i 


3-point 


64 Ge 


-4.80 


3.4 


-5.27 


3.2 


3-point 


68 Se 


-3.50 


3.0 


-4.55 


2.8 


3-point 


72 Kr 


-2.50 


2.6 


-4.68 


2.9 


3-point 


76 Sr 


-0.50 


1.2 


-4.33 


2.7 


3-point 


80 Zr 


-2.68 


2.8 


-4.61 


3.0 


5-point 


64 Ge 


-9.04 


4.3 


-8.31 


3.7 


5-point 


68 Se 


-9.95 


4.6 


-7.85 


3.6 


5-point 


72 Kr 


-9.78 


4.7 


-8.16 


3.7 


5-point 


7B Sr 


-4.53 


3.3 


-7.70 


3.6 


5-point 


80 Zr 


-7.75 


4.4 


-8.20 


3.9 



TABLE X: Same as in Table QX] for protons. 



Formula Nucleus 



HFBCS 



HTDA 







-E'cond 


Y U i v i 


E CO nd 


Y U i V i 


3-point 


64 Ge 


-3.55 


3.1 


-3.90 


2.9 


3-point 


68 Se 


-2.22 


2.5 


-3.15 


2.4 


3-point 


72 Kr 


-1.00 


1.7 


-3.40 


2.6 


3-point 


76 Sr 


-0.00 


0.0 


-2.94 


2.3 


3-point 


80 Zr 


-1.42 


2.1 


-3.01 


2.4 


5-point 


64 Ge 


-6.73 


3.9 


-6.12 


3.4 


5-point 


68 Se 


-7.34 


4.2 


-5.58 


3.2 


5-point 


72 Kr 


-7.46 


4.3 


-6.11 


3.4 


5-point 


76 Sr 


-2.36 


2.5 


-5.33 


3.1 


5-point 


80 Zr 


-5.70 


4.0 


-5.48 


3.3 



with a G- value fitted to the 3-point pairing indicator, no 
superfluid solution is found (see also Table IIIII for the 
values of pairing gaps). 



VII. PROTON-NEUTRON PAIRING IN THE 
HTDA APPROACH 

In spite of the wide recent theoretical interest paid 
to the proton-neutron pairing mode, it is still uncertain 
what is the exact importance of its T = component. 
Moreover, its connection with the Wigner energy is not 
completely clarified and there are a lot of controversies 
about other signatures. Since the T — pairing is ne- 
glected in all fits of the effective interactions in use to 
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TABLE XI: Ground-state properties of the five considered nuclei as functions of the x value (see text for details). 



Formula Nucleus 



r m (fin) Q20 (b) Q 4 o (b 2 ) 

a; = 0a; = lj: = 2 2: = 0a;=l x = 2 x = x = 1 2; = 2 

3-point 64 Ge 3.916 3.916 3.917 2.691 2.693 2.697 0.0065 0.0064 0.0062 

3-point 68 Se 4.012 4.012 4.012 -3.309 -3.308 -3.308 0.0523 0.0522 0.0522 

3-point 72 Kr 4.119 4.120 4.120 -4.601 -4.602 -4.603 0.1081 0.1082 0.1084 

3-point 76 Sr 4.239 4.240 4.240 7.499 7.501 7.506 0.2354 0.2356 0.2360 

3-point 80 Zr 4.320 4.320 4.321 8.446 8.446 8.450 0.1798 0.1795 0.1792 

5-point 64 Ge 3.917 3.918 3.918 2.688 2.690 2.696 0.0069 0.0067 0.0065 

5-point 68 Se 4.013 4.013 4.013 -3.308 -3.308 -3.307 0.0525 0.0524 0.0523 

5-point 72 Kr 4.120 4.120 4.120 -4.590 -4.592 -4.490 0.1073 0.1075 0.1078 

5-point 76 Sr 4.240 4.240 4.240 7.495 7.498 7.506 0.2351 0.2354 0.2361 

5-point 80 Zr 4.321 4.321 4.322 8.442 8.444 8.450 0.1811 0.1806 0.1801 



calculate the mean-field, this missing contribution may 
cause some artificial bias in the outcome. In the cal- 
culations of masses that make use of a macroscopic en- 
ergy (liquid-drop models), this problem is circumvented 
by adding an ad hoc Wigner term. It is however highly 
predictable that the inclusion of proton-neutron pairing 
correlations within a HFB approach would lead to novel 
features of the mean field, as for instance different de- 
formation properties since the usual spin-triplet (T = 0) 
pairing mode would tend to break the axial symmetry. 
However, since we do not carry out self-consistent HTDA 
calculations, the mean field itself is not affected by the 
presence of these proton-neutron correlations. 

Expectation values of various observables (such as the 
radii and quadrupole moments) specifying the correlated 
wave function have been calculated from solutions ob- 
tained after diagonalization of a full isoscalar and isovec- 
tor residual interaction. For the reasons mentioned in the 
previous paragraph, these quantities are not expected to 
change much in this perturbative treatment. In contrast, 
the change in energy brought in by the consideration of a 
full T z = part of the interaction is expected to be more 
significant. In the HTDA framework as applied here, and 
as far as relative variations are concerned, the correlation 
energy defined in Eq. (|2"5)) is relevant. 

The calculations reported in this section are performed 
in the 4p-4h space of pair excitations containing all new 
configurations that result from the coupling of neutron 
and proton states to produce the + ground state (includ- 
ing aligned proton- neutron pairs). The proton- neutron 
configurations considerably enlarge the dimension of the 
Hamiltonian matrices to be computed and diagonalized 
(up to ~ 10 5 ). As a result, the computing time becomes 
an issue. This is why we do not test here the conver- 
gence of the particle-hole expansion by going up to order 
6 as we have done in the previous section. A more de- 
tailed study of the proton- neutron pair correlations in the 



HTDA approach will be given in a forthcoming publica- 
tion. In addition, no evidence from previous theoretical 
approaches has been found for the T = collectivity. On 
the contrary, a rather strong quenching of the isoscalar 
pairing has been observed in the particular case of pf- 
shell nuclei. It is thus likely that we only have to deal 
with a somewhat weak neutron-proton pairing and con- 
sequently the 4p-4h limitation should not constitute a 
stringent constraint. 

Based on isospin invariance arguments we choose the 
strength V^ 1 = 1/2(V^ =1 + V^ =1 ) for the part of the 
residual interaction acting on neutron-proton two-body 
states. Since the actual T = pairing strength is un- 
known, in other words since one does not know what 
should be the data pertaining to the determination of 
a phenomenological such interaction, we adopt an ex- 
ploratory approach by varying the ratio x = V^f =0 /V^f =1 
of the residual interaction in the two isospin channels 
from 0.5 to 2.0 by steps of 0.5. 

In Table IXII the resulting GS deformations and radii 
are indicated only for three values of x, namely 0, 1 and 
2. As expected, these quantities do not vary significantly 
with x. 

The relative correlation energy with respect to the 
T = pairing mode, i.e., the difference between the val- 
ues of -Ecorr calculated with a given x value and with 
x = 0, is plotted in Fig. [3] as a function of x. The cor- 
relation energy induced by the T = mode is rather im- 
portant for large x values, so it is highly desirable to get 
some guidance on the actual value of x. One could argue 
that the shell- model estimates of Rcf. provide a value 
of about 1.5. However, the shell-model pairing definition 
with only (J,T) = (0,1) and (J,T) = (1,0) couplings 
is not consistent with the HTDA one where the S-iorce 
in use includes all multipolarities. Therefore, a system- 
atic comparative study of pair correlations in large shell- 
model calculations and in the HTDA approach would be 
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FIG. 3: Correlation energy (in MeV) normalized to the so- 
lution without T = pair correlations (x — 0) as a function 
of x for all considered nuclei. The case (a), respectively (b), 
corresponds to calculations with the pairing strength adjusted 
through the 3-point formula, respectively 5-point formula, in 
the T — 1 channel. 



very helpful to better assess a realistic value of x. 

Finally it is important to add that the energy shift due 
to the T = pair correlations depends on the pairing 
intensity in the T = 1 channel. Indeed these pair corre- 
lations are larger when the T = 1 pairing interaction is 
stronger. This suggests that, in a full treatment of these 
correlations, the adjustment of the pairing parameters in 
the T = 1 channel should be done with proton-neutron 
correlations. 



VIII. CONCLUSIONS AND PERSPECTIVES 

The purpose of this paper is to provide a firm basis 
to study the correlations present in deformed even-even 
N = Z nuclei in the A <~ 70 region. We mean, first 
of all, pairing correlations in the sense of Cooper pair 
excitation of nuclcons belonging to similar orbits. For 
these particular nuclei, proton-neutron correlations must 
clearly be included in the scheme. Moreover, RPA corre- 
lations should be further added. In view of the difficulties 
to be overcome in achieving this program, we have under- 



taken here initial steps by including only proton-neutron 
pairing correlations in an exploratory way. 

Since the various traditional pairing schemes often lead 
to a weak coupling regime, the Bogoliubov vacuum ansatz 
for the correlated ground state is rather inadequate. This 
is especially true for the not-so-heavy nuclei considered 
here. In contrast the HTDA framework represents a con- 
sistent and more physical way of handling various kinds 
of correlations, including RPA correlations, on the same 
footing. In the present work, we have chosen to use it 
from as good as possible HFBCS solutions with the usual 
pairing treatment which involves the T = 1 channel only. 

We have had, therefore, to carefully study the pair- 
ing channel which includes in particular the choice of a S 
interaction strength consistent with the data on atomic 
masses in the studied region. As a result, we have calcu- 
lated some deformation properties which are naturally of 
paramount importance to grasp the single-particle spec- 
troscopic properties. Upon determining the HTDA cor- 
related ground state, we have shown that the diffuseness 
of the Fermi surface is somewhat similar to what HFBCS 
yields but less sensitive to the fluctuations of the level 
density as a function of the mass number. 

In the last step, where the full T z = residual inter- 
action has been considered, we have established the fea- 
sibility of such calculations. Then, we have assessed in a 
quantitative way how much the relative amount of T = 
and T = 1 components of the residual interaction influ- 
ences the correlations properties. Even though in this 
study the T = 1 interaction strength is the one adjusted 
in the absence of a T z = component, it clearly appears 
from our results that a determination of the above ratio 
of isospin components is badly needed. 

The directions of improvement are therefore easy to 
perceive. First they imply a better understanding of 
the residual interaction. Then, RPA correlations should 
be included in a consistent way, which is currently in 
progress. 
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APPENDIX: SPECIFIC QUANTITIES USED TO 
ASSESS THE DEFORMATION PROPERTIES OF 
OUR SOLUTIONS 

In the HFBCS case, the expectation value of a local 
one-body operator may be expressed as a space integral 



12 



involving the local (diagonal in r) one-body reduced den- 
sity matrix. For the mass or isoscalar (neutron plus pro- 
ton) distribution, one defines the root-mean-square ra- 
dius r m and the quadrupole Q20 and hexadecapole Q40 
moments as 



with 



J cPr p(r) r 2 
A 



(31) 
(32) 



Q 2 o = J d?vp(v) (2z 2 -x 2 ~y 2 ) 1 

Q40 = / d 3 rp(r)r 4 r 4 °(0), (33) 

where Yg denotes the spherical harmonic of order I and 
magnetic quantum number 0. 

We then consider the equivalent spheroid which has the 
same root-mean-square radius and quadrupole moment 
as the actual nucleus. Denoting the semi-axes along the 
symmetry axis and perpendicular to it by a and c, re- 
spectively, we have 



Ar 2 = \{2a 2 + c 2 ), 
5 

O20 = \A{c 2 -a 2 ). 
5 



(34) 
(35) 



The fa parameter is then calculated for this equivalent 
spheroid by expanding the nuclear radius in polar coor- 
dinates according to the /^-parameterization [13| 



R(9) = 



y/l — a cos 2 9 

/ 00 N 



(36) 
(37) 



1-2- 

c 2 



(38) 



This allows us to calculate analytically the expression of 
02 for the equivalent spheroid as a function of a as 



2a 



_3_ 

2a 



Arcsin yfa. 



ae]0;l[ 
a = 



a < 



(39) 



1=1 



In the HTDA case, the above quantities have to be eval- 
uated in the correlated state and one cannot use any- 
more the usual generalized Wick theorem. Indeed, one 
has to evaluate matrix elements between two Slater de- 
terminants, generally different. One might keep, however, 
the Wick theorem framework by using mixed densities a 
la Lowdin [Lsl ] . Instead, here we reduce these many-body 
matrix elements into matrix elements evaluated between 
single-particle states, which makes the HTDA calculation 
of the above expectation values similar to the HFBCS 
ones. 
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